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Abstract. We show that growing neutrino models, in which the growing neutrino mass stops the 
dynamical evolution of a dark energy scalar field, lead to a substantial neutrino clustering on the 
scales of superclusters. Nonlinear neutrino lumps form at redshift z«l and could partially drag the 
clustering of dark matter. If observed, large scale non-linear structures could be an indication for a 
new attractive force stronger than gravity and mediated by the 'cosmon' dark energy scalar field. 
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Realistic quintessence models need to explain the recent increase of £1$ (z) by a factor 
of around ten or more [1, 2, 4, 3]. It has been recently proposed [6, 7, 8] that a growing 
mass of the neutrinos may play a key role in stopping the dynamical evolution of the dark 
energy scalar field, the cosmon. For a slow evolution of the cosmon, the scalar potential 
acts like a cosmological constant, such that the equation of state of dark energy is close 
to w = — 1 and the expansion of the universe accelerates. In these models, the onset 
of accelerated expansion is triggered by neutrinos becoming non relativistic. For late 
cosmology, z £ 5, the overall cosmology is very similar to the usual ACDM concordance 
model with a cosmological constant. 

An efficient stopping of the cosmon evolution by the relatively small energy density of 
neutrinos needs a cosmon-neutrino coupling that is somewhat larger than gravitational 
strength. This is similar to mass varying neutrino models [9, 10, 11, 12, 13, 14, 15, 16, 
17, 18, 19, 20, 21, 22, 29, 30], even though the coupling in those models is generically 
much larger. In turn, the enhanced attraction between neutrinos leads to an enhanced 
growth of neutrino fluctuations, once the neutrinos have become non-relativistic [10, 
12, 6, 23]. In view of the small present neutrino mass, m v (to) < 2.3eV, and the time 
dependence of m v , which makes the mass even smaller in the past, the time when 
neutrinos become non-relativistic is typically in the recent history of the universe, say 
zr ~ 5. Neutrinos have been free streaming for z > zr, with a correspondingly large free 
streaming length. Fluctuations on length scales larger than the free streaming length are 
still present at zr, and they start growing for z< zr with a large growth rate. This opens 
the possibility that neutrinos form nonlinear lumps [6, 24] on supercluster scales, thus 
opening a window for observable effects of the growing neutrino scenario. 

We show that neutrino perturbations indeed grow non-linear in these models. Non- 
linear neutrino structures form at redshift z ~ 1 on the scale of superclusters and beyond. 
One may assume that these structures later turn into bound neutrino lumps of the type 
discussed in [24]. Our investigation is limited, however, to linear perturbations. We can 
therefore provide a reliable estimate for the time when the first fluctuations become non- 
linear. For later times, it should only be used to give qualitative limits. 
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A crucial ingredient in this model is the dependence of the neutrino mass on the 
cosmon field 0, as encoded in the dimensionless cosmon-neutrino coupling j8(0) = 
~ dl d^ v ■ ^ or mcreasm g an d /3 < the neutrino mass increases with time m v = 

m v e~^^, where m v is a constant and j3 = j3 + <9j8/<9ln^>. The coupling j8(0) can 
be either a constant [6] or, in general, a function of 0, as proposed in [7] within a 
particle physics model. The cosmon field is normalized in units of the reduced Planck 
mass M = (8xGn)~ 1 / 2 , and j8 ~ 1 corresponds to a cosmon mediated interaction for 
neutrinos with gravitational strength. For a given cosmological model with a given time 
dependence of 0, one can determine the time dependence of the neutrino mass m v (t). 
For three degenerate neutrinos the present value of the neutrino mass m v (to) can be 
related to the energy fraction in neutrinos via £l v (to) = (3m v (fo)) / (94eVh ). 
We choose an exponential potential [1, 2, 4, 3]: 

y(0)=M 2 C/(0)=M 4 e - a ^ , (1) 

where the constant a is one of the free parameters of our model. 

The homogeneous energy density and pressure of the scalar field are defined in the 
usual way [8]. We can express the conservation equations for dark energy and growing 
matter in the form of neutrinos as [25, 26] 

p; = -3^(l+w^+ J 8(0)0'(l-3w v )pv , (2) 
p; = -3^(l+wv)Pv-/3(0)f(l-3wv)p v . 

The sum of the energy momentum tensors for neutrinos and the cosmon is conserved, 
but not the separate parts. We neglect a possible cosmon coupling to Cold Dark Matter 
(CDM), so that p' c = -3J?p c . 

For a given potential (1) the evolution equations for the different species can be 
numerically integrated, giving the background evolution shown in FIG.l (for constant 
j8) [6]. The initial pattern is a typical early dark energy model, since neutrinos are still 
relativistic and almost massless. Radiation dominates until matter radiation equality, 
when CDM takes over. Dark energy is still subdominant and falls into the attractor 
provided by the exponential potential (see [25, 26] for details). As the mass of the 
neutrinos increases with time, the term ~ /3p v in the evolution equation for the cosmon 
(2) starts to play a more significant role, kicking out of the attractor as soon as 
neutrinos become non-relativistic. This resembles the effect of the coupled dark matter 
component in [27]. Subsequently, small decaying oscillations characterize the — V 
coupled fluid and the two components reach almost constant values. The values of the 
energy densities today are in agreement with observations, once the precise crossing 
time for the end of the scaling solution has been fixed by an appropriate choice of the 
coupling /3 . At present the neutrinos are still subdominant with respect to CDM, though 
in the future they will take the lead (see [6] for details on the future attractor solution for 
constant /3). 

The evolution equations for linear perturbations (in Fourier space), in Newtonian 
gauge (in which the non diagonal metric perturbations are fixed to zero) [28], is fully 
described in [8]. 
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FIGURE 1. Energy densities of neutrinos (dashed), cold dark matter (solid), dark energy (dotted) and 
photons (long dashed) are plotted vs redshift. For all plots we take a constant j3 = —52, with a = 10 and 
large neutrino mass m v = 2.1 1 eV. 
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FIGURE 2. Longitudinal density perturbation for CDM (solid), v (dashed) and <j) (dot-dashed) vs 
redshift for k = O.lh/Mpc (upper panel) and k = l.lh/Mpc (lower panel, A = %Mpc). The neutrino 
equation of state (dotted) is also shown. The long dashed line is the reference ACDM. 



We numerically compute the linear density perturbations both using a modified ver- 
sion of CMBEASY [34] and, independently, a modified version of CAMB [37]. We plot 
the density fluctuations 5 ; as a function of redshift for a fixed k in FIG. 2. The neutrino 
equation of state is also shown, starting from 1/3 when neutrinos are relativistic and 
then decreasing to its present value when neutrinos become non relativistic. The turning 
point marks the time at which neutrino perturbations start to increase. At the scale of 
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k = QAh/Mpc (corresponding to superclusters scales) shown in FIG. 2a, neutrino per- 
turbations eventually overtake CDM perturbations and even force perturbations to in- 
crease as well, in analogy with dark energy clustering expected in [35, 36] within scalar 
tensor theories. Notice, however, that the scale at which neutrinos form nonlinear clumps 
depends on the model parameters, in particular the coupling /3, the potential parameter 
a and the present days neutrino mass. Those are related to the neutrino free-streaming 
length, the range of the cosmon field and its mass. A detailed investigation of the pa- 
rameter space will be performed in future work, but we mention that the model [7] with 
varying /3 gives qualitatively similar results. 




FIGURE 3. Redshift of first non linearities vs the wavenumber k for CDM (solid), v (dashed) and 
(dot-dashed). We also plot CDM for a reference ACDM model (dotted). 

We emphasize again that once the neutrinos form strong nonlinearities - neutrino 
lumps - one expects that nonlinear effects substantially slow down the increase of <5 V 
and even stop it. The magnitude of the CDM-dragging by neutrinos is therefore not 
shown - it might be much smaller than visible in FIG. 2. These remarks concern the 
quantitative interpretation of all the following figures, which are always computed in the 
linear approximation. 

Nevertheless, the linear approximation demonstrates well the mechanisms at work. 

In FIG. 3 we plot the redshift z„i at which CDM, neutrinos and § become nonlinear as a 
function of the wavenumber k. The case of CDM in the concordance ACDM model with 
the same present value of CI* and for massless neutrinos is also shown for reference 
(dotted line). The redshift z n i roughly measures when nonlinearities first appear by 
evaluating the time at which 8{z n i) = 1 for each species. The curves in FIG. 3 are 
obtained in the linear approximation, such that only the highest curves are quantitatively 
reliable. This concerns CDM for large k and neutrinos for small k. The subleading 
components are influenced by dragging effects and may be, in reality, substantially 
lower. 

We can identify four regimes: i) At very big scales (larger than superclusters) the 
universe is homogeneous and perturbations are still linear today, ii) The range of length 
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scales going from 14.5 Mpc to about 4.4 x 10 3 Mpc appears to be highly affected by the 
neutrino coupling in growing matter scenarios: neutrino perturbations are the first ones 
to go nonlinear and neutrinos seem to form clumps in which then both the scalar field and 
CDM could fall into. Note that the effect of the neutrino fluctuations on the gravitational 
potential induces CDM to cluster earlier with respect to the concordance ACDM model, 
where CDM is still linear at scales above ~87 Mpc. iii) For lengths included in the 
range between 0.9 Mpc and 14.5 Mpc, CDM takes over. That is in fact expected since 
neutrinos start to approach the free streaming scale. In this regime CDM drags neutrinos, 
and this effect may be overestimated in the linear approximation. Notice that in our 
model CDM clusters later than it would do in ACDM. There are two reasons for this 
effect. At early times, the presence of a homogeneous component of early dark energy, 
Clj, ~ 3/a 2 , implies that Q. m is somewhat smaller than one and therefore clustering is 
slower [38]. At later times, Q. m is smaller than in the ACDM model since for the same 
Q-tf,, part of 1 — £1$ = Q. m + £l v is now attributed to neutrinos. In consequence massive 
neutrinos reduce structure at smaller scales when they do not contribute to the clumping. 
The second effect is reduced for a smaller present day neutrino mass, iv) Finally, at very 
small scales (below clusters), CDM becomes highly non linear and neutrinos enter the 
free streaming regime, their perturbations do not growth and remain inside the linear 
regime. 

Maximum neutrino clustering occurs on supercluster scales and one may ask about 
observable consequences. First of all, the neutrino clusters could have an imprint on 
the CMB-fluctuations. Taking the linear approximation at face value, the ISW-effect of 
the particular model presented here would be huge and strongly ruled out by observa- 
tions. However, non-linear effects will substantially reduce the neutrino-generated grav- 
itational potential and the ISW-effect. Further reduction is expected for smaller values 
of /3 (accompanied by smaller a). It is well conceivable that realistic models for the 
growing neutrino scenario lead to an ISW-effect in a range interesting for observations. 

A second possibility concerns the detection of nonlinear structures at very large length 
scales. Such structures can be found via their gravitational potential, independently of 
the question if neutrinos or CDM source the gravitational field. Very large nonlinear 
structures are extremely unlikely in the ACDM concordance model. An establishment 
of a population of such structures, and their possible direct correlation with the CMB- 
map [39, 40, 41, 42, 43], could therefore give a clear hint for "cosmological actors" 
beyond the ACDM model. For any flat primordial spectrum the gravitational force will 
be insufficient to produce large scale clumping, which could thus be an indication for a 
new attractive force stronger than gravity - in our model mediated by the cosmon. 
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